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ABSTRACT 

We use semi- analytical modelling of galaxy formation to predict the redshift-size-evolution of el- 
liptical galaxies. Using a simple model in which relative sizes of elliptical galaxies of a given mass 
correlate with the fraction of stars formed in a star burst during a major merger event, we are able 
to reproduce the observed redshift-size-evolution. The size evolution is a result of the amount of cold 
gas available during the major merger. Mergers at high redshifts are gas-rich and produce ellipticals 
with smaller sizes. In particular we find a power-law relation between the sizes at different redshifts, 
with the power-law index giving a measure of the relative amount of dissipation during the mergers 
that lead to the formation of an elliptical. The size evolution is found to be stronger for more massive 
galaxies as they involve more gas at high redshifts when they form, compared to less massive ellipti- 
cals. Local ellipticals more massive than 5 x 10 11 M© will be approximately 4 times larger than their 
counterparts at z = 2. Our results indicate that the scatter in the size of similar massive present day 
elliptical galaxies is a result of their formation epoch, with smaller ellipticals being formed earlier. 
Subject headings: galaxies: elliptical - galaxies: interaction- galaxies: structure - galaxies: evolution 
- methods: numerical 



1. INTRODUCTION 

Elliptical galaxies have been shown to have formed the 
bulk of their s tars at high redshifts within an intensiv e 
star burst (e.g. lCimatti et alJl2004 iThomas et"afl [2005L 
supporting a scenario in which elliptical galaxies formed 
in a 'monolithic' collapse and then continued evolving 
passively. On the other hand, many stellar dynami- 
cal properties of elliptical galaxies can be explained by 
the merger of similar massive galaxies (e.g. Barnes & 
Hernquist 1992; Naab & Burkert 20 03; Naab, Jesseit fc 
Burke rt 2006) as initially proposed bv lToomre fc Toomrd 
(fl972(h The predicted merger rate of galaxies within the 
CDM paradigm and the observa tions are in fair agree- 
ment l)Khochfar fc Burkertl l2001L As to the mainly old 
stellar populations of elliptical galaxies, iDe Lucia et alJ 
(2006) show that these can be recovered within the 
merger scenario which comes as a consequence of massive 
elliptical galax ies forming by dry mergers between ellip- 
tical galaxies llKhochfar fc Burkerti r2003). In addition, 
it has been shown that dry mergers can indeed explain 
the kinematical properties of massive elliptical galaxies 
(Khochfar & Burkert 2005; Naab, Khochfar & Burkert 
2006). and might be able to even account for the par- 
tially depleted cores of luminous ellipticals (e.g. lGrahaml 
2004). 

The size evolution of elliptical galaxies could serve as 
an additional possible test for these two models. Re- 
cent size measurements of ellipticals carried out at high 
rcdshift indicate that they are mu c h sma l ler than their 
local counter parts ijDaddi et alJ 120051 iTruiillo et alJ 
120051) . Numerical simulations bv lNaab fc Truiillol l|2005F 
which use self-consistent cos mological orbital parameters 
ijKhochfar fc Burkerd l200rj|) to set up mergers between 
disk galaxies, find that the size of the remnant is of the 
same ord er as the size of the pro genitor disk. The simu- 
lations of lNaab fc Truiillol l)2005|) did not include any gas 
or star formation in contrast to simulations carried out 
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bv I Springe! & Hernquist (2005). The latter authors find 
that those stars which existed in the progenitor disks be- 
fore the merger, which we call the quiescent component, 
have a ~ 5.7 times larger effective radius than the stars 
which formed during the merger in a violent star burst, 
which we will call the merger component of the elliptical. 
However, the ratio of the effective radii of the two com- 
ponents is very likely to depend on the structure of the 
progenitors and on the way the merger takes place. The 
total effective radius of the remnant will depend on the 
mass fraction of each component, with a larger merger 
component leading to a smal l er rem nant. 

Recently, iKhochfar feljilkl ((20051 hereafter KS) inves- 
tigated the fraction of merger and quiescent components 
in early-type galaxies, finding that the quiescent compo- 
nent is decreasing with redshift and increasing with mass 
up to a characteristic mass scale Mc = 3 x 10 10 Af© 
(|Kauffmann et all 12003. hereafter K03) at which it be- 
comes constant. In this letter we follow up on their study 
and predict the size-evolution of elliptical galaxies within 
the CDM-paradigm. 

2. MODEL 

We use semi-analytical modelling of galaxy formation 
to predict the merger and quiescent components of el- 
liptical galaxies. The dark matter history is calculated 
using the merger tree proposed by iSomerville fc KolattJ 
(1999) with a mass resolution of 2 x 10 9 Mq . The baryonic 
physics within these dark matter halos is calcul ated fol- 
lowing recipes presented in lSpringel et alJ i2001l) includ- 
ing a model for the reionizing background bv lSomervilkj 
(2002). In our simulation, we assume that elliptical 
galaxies form whenever a major merger (M1/M2 < 3.5 
with Mi > Mi) takes place. We assume that during this 
process all the cold gas which was in the progenitor disks 
will be consumed in a central starburst, adding to the 
spheroid mass, and that all stars in the progenitor disks 
will be scattered into the spheroid too. Furthermore we 
allow the stars of satellite galaxies in minor mergers to 
also contribute to the spheroid. During the evolution of 



2 



n.o < ion m, < ii.5:i 



10.0 < log M. < 10.5 : 




1 1.5 < log M, < 12.0 



10.5 < log M. < 11.0 



0.5 




5 < log M, < 10.0 



0.5 



Fig. 1. — Histograms show the conditional probability density 
p(l — M q /Mi, u i\M*) of the merger component in spheroids as a 
function of stellar mass M, in six different mass bins. Solid lines 
show log-normal distributions fitted to the data. Results are shown 
for elliptical galaxies only. 

a galaxy, we keep track of the origins of all stars brought 
into the spheroid and attribute them to two categories, 
merger and quiescent, where the first incorporates stars 
formed during a starburst in a major merger and the lat- 
ter includes stars previously formed in a disk and added 
to the spheroid during a major merger. Each star will 
carry along its label and not change it, which means that 
if a star was made in a merger of two progenitor galaxies 
and the remnant of that merger participated in another 
merger, the star will still contribute to the merger com- 
ponent of the final remnant. 

For more modelling details, we refer the reader to KS 
and references therein. Throughout this paper, we use 
the following set of cosmological parameters: f^o = 0.3, 
n A = 0.7, n b /n Q = 0.15, <7 8 = 0.9 and h = 0.65. We 
note here that running our code with the la test cosmolog- 
ical p arameters from the WMAP mission l)Spergel et al.l 
2006) only changes our results slightly and that most of 
the effects are compensated for by the free model param- 
eters for the star formation efficiency and the supernova 
feedback. 

3. RELATIVE SIZES OF ELLIPTICALS 

As we have mentioned above, the relative effective radii 
of merger and quiescent components are very likely to be 
dependent on the physical properties under which the 
merger takes place. To minimise this effect, we compare 
the relative sizes for ellipticals of similar mass. By doing 
so, we can assume that e.g. the potential depth is the 
same and that feedback effects have the same efficiency. 
Furthermore, it has been shown that the orbital param- 
eters of merging halos are on ave rage independent of 
the re mnant halo mass and redshift ijKhochfar fc BurkertJ 
2006), which allows us to assume that merging will take 
place on average with the same orbital set-up. 

The SDSS study revealed that the size distribution of 
galaxies with the same mass is log-normal and that the 
varian ce o"i n fl e of the distribution is a function of mass 
fK03: IShen etalll200l hereafter S03). Above - 10 10 
Mq, the variance drops until it becomes approximately 
constant for galaxies more massive than ~ 10 11 Mq. S03 




Fig. 2. — Dependence of the variance a g ^ of the conditional 
probability density p(l — Mq/My u i\M<t*) as a function of stellar 
mass M„ of elliptical galaxies. The dotted line shows the fit using 
the empirical formula of S03. 

show that their model, in which elliptical galaxies form 
by continued merging of galaxies, could reproduce the 
scatter in the size distribution of massive elliptical galax- 
ies, by assuming the size distribution of the progenitors 
to be log-normal. However, the origin of the scatter in 
the size distribution of the progenitors remains not well 
understood. We start off by investigating the scatter in 
the merger component and comparing it to the scatter 
in the observed size distribution of early-type galaxies. 

In Fig. ^ wc show the distribution of the merger com- 
ponents found in elliptical galaxies of various masses M* , 
where elliptical galaxies are those having more than 65% 
of their stellar mass in their bulge component. We have 
fitted the data with log-normal distributions following 
K03 and S03 as: 
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with a, b and ov^ as free parameters, x = 1 — M q /Mb u i as 
the merger component in a spheroid of mass Mf, u i, and 
M q as the quiescent component of that spheroid. As can 
be seen, the log-normal distribution provides excellent 
fits to our simulated data. We can now try to compare 
our variance a s ^ as a function of galaxy mass to the ob- 
served variance <Ti n i? e in the SDSS. S03 give an empirical 
fitting formulae to their results, 



<r sb = cr 2 



<J\ — a 2 



i + (m*/m ) 2 



(2) 



with o"i, (72 and Mq as free parameters. The fit of Eq. 
12 to our data is shown in Fig. [21 Our data seems to be 
fitted well by Eq. |2Jand the trend in the observations can 
be recovered. The mass scale Mq, characteristic for the 
transition point from large scatter to small scatter in the 
distributions is best fitted by a value of Mq w 1.7x 10 10 
Mq in our simulations, which is a about a factor 2 smaller 
than the value of 3.89 x 10 10 M Q suggested by S03 but 
still in reasonable agreement. For completeness we also 
give the values of the other fitting parameters which are 
(7i = 0.31 and a 2 = 0.14. 
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Fig. 3. — The evolution of sizes for early-type galaxies with 
respect to the sizes of their local counter parts at a redshift z = 0.1. 
The upper panel shows the size-evolution for early-type galaxies 
larger than 3 X 10 10 fey 2 Mq and the lower panel for early-type 
galaxies larger than 6.6 X lO 10 /^ 2 Mq. 

The results presented here suggest that the scatter in 
merger components 1 — M q /Mb u i behaves like the scatter 
in the size distribution of elliptical galaxies. We now 
make the simplified assumption that for each mass bin 

hxil-MjM^ochiiRe), (3) 

where Re is in units of kpc and the proportionality con- 
stant reflects the average physical conditions that led to 
the formation of the elliptical galaxy and. We note that 
other effects may influence the relative sizes of elliptical 
galaxies too and that we here only focus on the contribu- 
tion to it by the merger component. Using the assump- 
tion of proportionality one can now calculate the relative 
sizes of elliptical galaxies of the same mass by knowing 
their different merger components and using: 

Re(zx) = Re /d {z ) with z x > z , (4) 

with d ee ln(l - M qfl /M bu i)/\n(l - M qA /M bul ) as the 
dissipation factor which gives a measure of the relative 
amount of dissipation that led to the formation of an 
elliptical. The scatter in the size distribution of ellipticals 
decreases with mass and later becomes constant. The 
reason for the same behaviour in the scatter of the merger 
component is that most massive ellipticals have their last 
major merger in a small redshift window not too far back 
in time. As a consequence the conditions regarding the 
gas fraction involved in the merger are very similar and 
the scatter is small. 

4. SIZE EVOLUTION OF ELLIPTICALS 

We now can test size evolution as a function of red- 
shift predicted by Eq. 0] for elliptical galaxies of a given 
mass and compare it to the observations. Since we can 
only predict relative sizes between elliptical galaxies of 
approximately the same mass, we will normalise sizes 
to the SDSS sample. We calcul ate the size evolutio n 
for the same redshifts presented in lTruiillo et alJ |2005). 
The authors took the mean effective radii of the ln(i?. c ) 
distribution for galaxies above two mass thresholds of 
3 x 10 lc V n 2 Mo and 6.6 x \Q w h^ M n from the SDDS 



Fig. 4. — The predicted evolution of sizes for early-type galaxies 
with respect to the sizes of their local counter parts at a redshift 
z = 0.1 divided into four different mass bins. 

sample of early-type galaxies and divided the effective 
radii of early-type galaxies at higher redshifts by this 
value. After arranging their galaxies in various redshift 
bins they calculated the means of these ratios and pre- 
sented these values. We here use the same method to 
compare our results to theirs. Our zero-point for individ- 
ual galaxies is taken to be the mean value of 1 — M q /Mb u i 
in our ln(l — M q /M}, u i) distribution for galaxies of the 
same mass at a redshift of z = 0.1. In Fig. [3] we show 
the expected evolution of sizes. For both cases of limit- 
ing masses, the agreement is excellent. It appears that 
the difference in sizes is more significant for massive early 
type galaxies. In Fig. 0we predict the size-evolution in 
four different mass ranges based on the relative amount 
of their merger component. While local early- type galax- 
ies between 10 10 M Q and 10 11 M are around 1.25 times 
larger than their counterparts at z = 2, local early types 
with masses larger than 5 x 10 11 M will be approxi- 
mately 4 times larger than their counterparts at z = 2. 
This dramatic change in sizes in our model results from 
massive galaxies at high redshifts forming in gas-rich 
mergers while galaxies of the same mass at low redshifts 
form from gas-poor mergers (KS). This size-evolution 
might be an overestimate as the modelled galaxies suffer 
from over cooling of gas (KS) which is likely to overes- 
timates the merger component at high redshift and the 
quiescent component at low redshift due to the shorter 
time between consecutive major mergers at high redshifts 
compared to low redshifts. 

5. DISCUSSION AND CONCLUSIONS 

We have connected the scatter in the merger compo- 
nents 1 — M q /Mbui with the scatter in the sizes of ellipti- 
cal galaxies. Following KS, the scatter in the merger 
components of elliptical galaxies is a result of differ- 
ent formation epochs. Ellipticals forming early have 
larger merger components, as they were formed in gas- 
rich mergers, hence ellipticals with smaller effective radii 
must have formed earlier. We tested this assumption 
by predicting the size-evolution of elliptical galaxies at 
different redshifts a nd the agreement with the data of 
Truiill o et al.l l)2005|) is excellent. It is important to note 
that we normalise the proportionality in our relation be- 
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tween the merger component and the effective radius by 
the observed scatter in the local galaxy sample, and then 
go ahead and predict how sizes at earlier redshifts com- 
pare to local sizes. 

Our results demonstrate that the strongest size evolu- 
tion is for massive elliptical galaxies. Local early-types 
with masses larger than 5 x 10 11 M Q will be approxi- 
mately 4 times larger than their counterparts at z = 2. 
This extreme size-evolution is a reflection of the progen- 
itors having larger gas fractions at high redshifts and 
becoming more 'dry' at low redshifts. In addition the 
low redshift progenitors have stellar disks that are more 
massive than the available amount of cold gas, hence in- 
creasing the quiescent fraction of the remnant. The pro- 
genitors however, are bulge dominated with a disk com- 
ponent not more massive than ~ 20% of the total mass, 
thus the size-evolution appears to be connected to the 
occurrence of bulge-dominated dry mergers with time. 
The most massive elliptical galaxies in our simulations 
undergo on average between one and two substantial dry 
mergers between z = 2 and today. 

We find that the relative amount of dissipation in- 
volved in the mergers relates to the size by a power-law 
as described in Eq. 0] where the power-law exponent d is 
the dissipation factor. If this relation holds one can try to 
measure the relative amount of dissipation by measuring 
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the relative sizes of recently formed ellipticals of the same 
mass at different redshifts. Even though theoretically it 
sounds straightforward to measure the dissipation fac- 
tor we acknowledge that it is observationally not an easy 
task. The main problem here will be to identify elliptical 
galaxies that just formed. This is important as e.g. con- 
tinued accretion of satellites or cold gas and subsequent 
star formation will alter the size of the merger remnant . 
One way of identifying recently formed ellipticals at high 
redshift might be by looking for signs of recent star for- 
mation. If there is sufficient gas involved in the major 
mer ger, simulations show that a starburst will be ignited 
fe.g iBarnes fc Hernauistlll992D whose signature may be 
measurable. 

Our results presented here support a picture in which 
ellipticals form in mergers. Future observations of high 
redshift ellipticals will allow one to make more accurate 
comparisons to the model we introduced here and will 
allow us to estimate the role of dissipational processes 
during major mergers of galaxies. 
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paper. 
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